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1 . INTRODUCTION 


JUr  f 

FRICKDIto  VAOkfo UKrNOT  riLMKD 


This  paper  started  as  a tutorial  exercise  to  help  tighten  up  some  of 
the  concepts  that  a number  of  groups  have  been  referring  to  somewhat 
loosely.  The  idea  was  to  start  with  a quite  general  statement  of  the  fire 
control  problem  as  in  equation  2.1  and  to  go  from  there  to  define  the  roles 
of  dither,  noise,  and  weapon  precision,  to  define  the  concept  of  optimal 
algorithms,  and  to  discuss  the  idea  of  optimal  maneuvers. 


All  that  is  contained  here,  but  in  addition  there  are  some  more  quan- 
titative results  as  in  equation  6.12  which  developed  as  the  paper  was 
being  written.  These  allow  bounds  to  be  put  on  the  ability  of  the  target 
to  minimize  the  best  that  the  fire  control  can  do,  and  they  allow  this 
characterization  to  be  expressed  by  quantities  calculated  from  only  the 
correlation  function  of  the  process  — whatever  the  process  --  Gaussian 
or  otherwise. 


2.  STATEMENT  OF  THE  PROBLEM 

The  mathematical  problem  analyzed  in  this  paper  is  posed  in  terms 
that  relate  it  to  fire  control  problems  for  gun  systems.  The  gunner  has 
data  on  the  behavior  of  the  target  over  some  subset  of  times  taken  from 
the  set  of  times  less  than  or  equal  to  tQ.  These  data  are  most  likely 
corrupted  by  noise.  At  time  tQ  the  gunner  fires.  The  bullet  is  committed 
to  its  trajectory  at  time  tQ  and  arrives  at  the  target  plane  at  time  t0 + T. 
The  target  plane  can  be  taken  as  the  plane  containing  the  bullet  and  a 
reference  point  on  the  target  and  perpendicular  to  the  relative  velocity 
vector  of  the  bullet  and  target. 

It  is  usually  convenient  to  consider  the  distribution  of  bullet  posi- 
tions in  the  target  plane,  and  this  gives  the  problem  an  essentially  two 
spatial  dimensional  character  — at  least  if  the  definition  of  the  target 
plane  can  be  made  unambiguously.  Since  the  prime  motivation  for  this 
paper  was  the  air  defense  gun  problem,  we  shall  emphasize  the  two-dimensional 
problem  --  when  the  dimensionality  is  important  to  the  discussion.  However, 
we  should  point  out  that  there  is  an  important  class  of  problems  that  are 
fundamentally  one-dimensional  — for  example,  the  case  in  which  a vehicle 
is  constrained  to  move  on  the  earth  and  that  vehicle  is  being  engaged  by  a 
gun  whose  trajectory  is  virtually  parallel  to  the  earth. 

The  gunner  may  well  fire  more  than  one  round  at  the  target.  In  our 
analysis  we  assume  that  he  does  not  use  data  on  previous  hits  and  misses 
to  correct  his  fire  — each  shot  is  handled  as  an  independent  operation  on 
the  part  of  the  gunner. 

Likewise,  the  pilot  of  the  target  performs  maneuvers  that  are  con- 
strained by  physical  restrictions  on  his  vehicle,  requirements  of  his 
mission,  and  his  desire  to  avoid  being  hit.  The  pilot  may  institute  a 
series  of  evasive  maneuvers  when  he  observes  the  firing  of  the  gun,  but 
he  does  not  cue  the  individual  twists  and  turns  of  his  maneuver  on  the 
firing  of  each  shot  from  the  gun. 
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We  assume  that  the  gunner  has  complete  knowledge  of  the  statistics 
of  the  target's  maneuvers.  We  also  assume  that  the  pilot  is  aware  of 
this  fact  and  also  knows  that  the  gunner  is  smart  enough  to  know  how  to 
optimize  his  fire  control  on  the  basis  of  this  statistical  knowledge. 

In  symbols,  the  problem  becomes  one  of  minimizing  p (by  the  pilot) 
and  maximizing  p (by  the  gunner)  where 

P = / b[q  - c - f (z)]  j (c)v(q -x)g(x,z)dqdcdxdm(z),  2.1 

and  where 

p 3 the  probability  of  killing  the  target  given  a shot  at 
time  tQ 

b(*)  ■ the  probability  density  function  of  the  ballistic 
dispersion  of  the  gun  about  the  aim  point 

q = the  vector  position  of  the  bullet  in  the  target  plane 

c = that  extra  dispersion  called  dither  that  is  sometimes 
added  to  the  aiming  of  the  gun 

f(*)  3 the  fire  control  algorithm.  This  is  a functional 
relating  the  past  performance  of  the  target  to  a 
prediction  of  the  future  position  of  the  target. 

z = y + n 

y = history  of  the  path  of  the  target  taken  over  some  sub- 
set of  the  set  of  times  less  than  or  equal  to  tQ. 

For  example,  this  may  be  the  history  of  the  target 
over  the  interval  t0  < t < t . Or  the  history  could 
be  taken  over  several  disjoint  intervals  as  the 
target  pops  in  and  out  of  view. 

n * the  noise  that  represents  errors  in  measuring  y.  This 
is,  of  course,  taken  over  the  same  subset  of  times 
as  is  y. 

j(»)  3 the  probability  density  function  of  dither 

v ( • ) 3 the  probability  that  a hit  at  the  point  described  by 
its  argument  will  produce  a kill.  Note  that  this 
is  not  a probability  density  function, 

x 3 the  future  position  of  the  target;  that  is,  the  posi- 
tion at  time  t0  + T 

g(*,*)  3 the  conditional  probability  density  function  of  x 
given  z 

m(z)  3 the  probability  measure  over  the  ensemble  Z 

Z 3 the  ensemble  of  functions  z 
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We  will  also  have  use  for  the  following: 

m(y)  ■ the  probability  measure  over  the  ensemble  Y 
Y * the  ensemble  of  possible  target  path  histories  y 
m(n)  » the  probability  measure  over  the  ensemble  N 
N * the  ensemble  of  noise  histories  n 

It  is  more  convenient  to  use  g(x,z)dx  than  the  more  complete  form 
dxG(x,z),  and  we  do  this  also  since  later  results  indicate  that  it  is 
important  to  the  pilot  to  minimize  sup  g(x,z). 

x 

Usually  the  extent  of  the  target  as  represented  by  v(«)  is  small 
compared  to  errors  in  predicting  the  future  position  of  the  target  — at 
least  in  those  cases  in  which  those  errors  are  large  enough  to  justify 
the  efforts  of  this  paper.  We  thus  replace  v(«)  by  A 6(*)  where  A is 
the  so-called  vulnerable  area. 

This  gives  us 

p = A J b[x  - c - f (z)]  j (c)g(x,z)dcdxdm(z) . 2.2 

We  shall  use  equation  2.2  in  all  that  follows. 

3.  THE  OPTIMAL  ALGOR ITttl 

If  everything  except  f(-)  is  given  in  equation  2.2,  then  the  gunner 
need  only  solve  the  problem  of  minimizing  p by  choosing  f(*)»  and  we  can 
write  equation  2,2  as 


p = J*  h[f  (z)  ,z]dm(z)  3.1 

by  formally  integrating  with  respect  to  all  the  other  dummy  variables. 

We  define  a set  of  functionals  F which  are  those  f(»)  such  that 
f (z)  is  measurable  in  Z.  Let  Z*  be  a subset  of  Z such  that  the  measure 
of  Z*  is  unity.  Equality  over  Z*  is  the  usual  equality  "almost  everywhere." 

We  assume  that  for  each  element  z £ Z*  there  is  an  x called  f(z) 
such  that 

h[$(z),z]  = sup  h(x,z)  3.2 

x 

and  we  assume  that  f(«)  C F.  Then 

p = J h[f  (z),z]dm(z)  > /h[f  (z),z]dm(z)  3.3 

for  all  f (• ) C F. 
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The  functional  f (•)  is  the  (perhaps  not  unique)  optimal  fire  control 
algorithm.  By  using  it,  the  gunner  takes  full  advantage  of  the  informa- 
tion available  to  him.  He  can  do  no  better.  The  remainder  of  his  error 
is  in  the  category  of  completely  unpredictable. 

Throughout  this  paper  we  shall  assume  that  the  gunner  should  maxi- 
mize the  individual  kill  probability  for  each  round  he  fires.  This 
assumption  is  not  strictly  true  for  rounds  fired  in  "bursts."  It  is  a 
satisfactory  assumption  for  rounds  fired  at  times  that  are  sufficiently 
different  that  the  target's  locations  at  the  two  times  are  not  strongly 
correlated.  It  is  also  valid  if  the  individual  rounds  have  a low  proba- 
bility of  killing  the  target  given  a hit.  However,  the  problem  with 
the  assumption  is  not  associated  with  predicting  the  target's  future 
position  but  rather  is  one  of  "overkilling"  the  target.  That  is,  if 
the  gunner  has  a very  effective  round,  he  would  prefer  to  fire  a second 
round  at  a somewhat  different  place  than  his  first  round  — if  the  first 
round  hit  the  target,  the  target  was  killed;  if  the  first  round  missed, 
why  send  a second  round  there? 

We  shall  sav  more  about  this  issue  after  we  have  some  of  the  argu- 
ments in  the  ne  section  and  have  established  the  idea  of  the  relative 
merit  of  deterrr  istic  algorithms.  Nevertheless,  we  must  point  out 
that  for  gunfire  against  aircraft  we  are  perhaps  concerned  about  not 
overlapping  a few  square  meters  of  vulnerable  area  while  irreducible 
errors  in  prediction  may  spread  the  target  over  hundreds  of  square 
meters.  Also  we  should  point  out  that  estimating  the  kill  of  a burst 
of  I rounds  as 

Pj  - 1 - U-P)1  3.4 

is  an  overestimate  and  as  such  is  consistent  with  the  strategy  of  our 
analysis  in  this  paper. 

4.  THE  ROLES  OF  DITHER  AND  GUN  ACCURACY 

In  this  section  and  in  the  next,  we  shall  "peel  the  onion"  from 
the  point  of  view  of  the  gunner.  That  is,  we  shall  simplify  equation  2.2 
by  optimizing  those  things  over  which  the  gunner  has  control  or  over 
which  he  could  have  control  if  he  were  to  choose  to  improve  his  weapon 
system. 

As  mentioned  above,  we  shall  do  this  analysis  with  the  assumption 
that  the  optimal  strategy  for  the  gunner  is  to  maximize  p for  each 
round.  This  is  probably  an  adequate  assumption  for  our  purposes,  but  for 
completeness  at  least  we  shall  then  show  that  the  ideas  of  this  section 
can  be  generalized  if  this  assumption  is  inadequate. 

First  consider  dither.  We  shall  show  that  dither  is  less  effective 
than  using  the  actual  target  statistics.  To  do  this,  we  let 


p = A j b[x  - f (z)]g(x,z)dxdm(z),  4.1 

where  f(«)  is  optimal  for  the  dither less  case.  That  is,  let 

A J*b[x  - f (z)]g(x,z)dxdm(z)>  A J b[x  - f (z)]g  (x, z)dxdm  (z)  4.2 

for  all  f € F. 

Then  if  we  let  f(z)  + c = fc(z)  and  note  that  fc(z)  £ F,  we  have  by 
integrating  both  sides  of  equation  4.2  multiplied  by  j (c) 

- c - f (z)]g(x,z)dxdm(z)  | j (c)  dc 

- J"tp^  j (c)dc  = p 4.3 


since  p is  independent  of  c and  j(*)  is  a probability  density  function. 

Thus  "the  gunner  does  at  least  as  well  by  using  the  optimal  predic- 
tion algorithm  in  the  ditherless  case  as  he  does  by  using  dither  and  the 
same  or  any  other  prediction  algorithm."  Dither  has  its  role,  however, 
if  the  target's  statistics  are  unknown  or  for  some  other  reason  a less 
than  optimal  algorithm  is  used.  This  point  is  illustrated  in  Appendix  A. 

For  our  main  effort,  however,  we  shall  assume  that  our  "smart"  gunner 
knows  g(*,*)  and  uses  j (• ) = 60)  and  optimizes 


p = A f b[x  - f (z)]g(x,z)dxdm(z) . 


We  can  argue  similarly  that  accurate  guns  are  better  than  less 
accurate  guns.  In  particular,  if  b(*)  is  the  convolution  of  bjO)  and 
b20)»  that  is 


and  if 


b(C)  = / bj  (£  - n)b2 (n)dn. 

Pi  = A / bjjx  - fi(z)]gCx,z)dxdm(z) 


is  optimized  for  each  i = 1 and  2,  then 


P < Pi 


for  each  i where  p is  optimal  for 


p = A / b[x  - f (z)]g(x,z)dxdm(z) . 
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The  proof  follows  the  lines  of  the  proof  on  dither  by  substituting 
equation  4.5  into  4.3  and  treating  one  of  the  bjf*)  as  b(*)  and  the  other 
as  j(«)  in  equations  4.1,  4.2,  and  4.3. 

Two  special  cases  are  important.  First,  since 
b(£)  = / b(£  - n)6(n)dn, 

the  gunner  does  best  of  all  with  a perfect  gun.  So  we  will  take  another 
layer  from  the  onion  and  henceforth  we  shall  use  b(*)  = 6(*)  as  the  best 
course  of  action  for  the  gunner.  This,  by  the  way,  does  not  represent  a 
physically  impossible  situation.  Guns  that  are  well  designed  have  accu- 
racies that  are  usually  small  compared  with  the  other  errors  of  the  weapon 
system  — so  much  so  in  fact  that  additional  dispersion  (dither)  is  often 
added  to  eliminate  biases. 

Now  we  want  to  optimize  (or  rather  the  gunner  wishes  to  optimize) 
p = A /g[f  (z)  ,z]dm(z)  4.9 

A 

which  he  does  by  picking  f(«)  so  that 

g[f(z),z]  > g[f(z),z]  4.10 

for  all  f(-)  £ F and  almost  everywhere  in  z.  Then  we  have 

P = A J g[f  (z)  ,z]dm(z)  . 4.11 

Second,  since  Gaussian  functions  are  usually  good  approximations  for 
ballistic  dispersion,  and  since  if  b(-),  bj(.),  b2(.)  are  Gaussian  with 
standard  deviations  a,  Oj , and  a2,  and 

a2  = Oj2  + ct22  4.12 

then 

b(£)  = f bj (5  - n)b2(n)  dn,  4.13 

it  follows  that  a lower  dispersion  gun  (c^  or  a2)  is  better  than  a higher 
dispersion  gun  (a)  if  the  optimal  prediction  algorithm  for  b^*)  or 
b2(0  is  used  for  the  lower  dispersion  gun  and  any  algorithm  is  used  for 
the  higher  dispersion  gun. 

Now  consider  a burst  of  rounds  where  each  round  is  highly  lethal 
(a  hit  is  a kill)  and  such  that  the  burst  is  fired  before  the  first  round 
reaches  the  target.  Certainly  the  gunner  should  fire  the  first  round  at 
x0  = f (z)  as  in  equation  4.10.  We  can  show  by  arguments  similar  to  that 
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above  that  a second  round  fired  at  time  tj  should  be  aimed  at  the  position 
Xi  which  maximizes  g(x,xo,Zj),  the  conditional  distribution  of  the  target 
at  time  tj  + T given  that  the  target  was  not  at  xQ  (within  the  size  of 
the  target)  at  time  tQ  + T and  given  the  history  Zj  up  to  time  tj. 

Likewise  a third  round  fired  at  time  t2  should  be  aimed  at  the  posi- 
tion which  maximizes  g (x,Xj ,x0, z2) , the  conditional  distribution  of  posi- 
tions of  the  target  at  time  t2  + T given  that  the  target  was  not  at  Xj  at 
time  tj  + T and  was  not  at  x0  at  time  t0  + T and  given  the  history  z2  up 
to  time  t2.  Thus  we  can  define  an  optimal  algorithm  which  is  determinis- 
tic. That  may  be  of  more  theoretical  than  practical  interest  — but  there 
it  is. 


5.  NOISE 

As  for  dither,  noise  is  degrading  if  the  optimal  prediction  algorithm 
is  used  for  the  noiseless  case.  To  show  this,  we  rewrite  equation  4.9  in 
the  form 


P = A / gt [f (y  ♦ n) ,y]dm(y)dm (n) , 5.1 

where  gjC*,’)  is  the  conditional  probability  density  function  of  the 
future  position  given  the  actual  history  of  the  target' s^path.  The  argu- 
ment of  f(-)  is,  of  course,  the  noisy  history.  Now  let  f(-)  be  such  that 

gJffyD.y]  > gJfM.y]  5.2 

almost  everywhere  in  y.  Now,  letting  fn(y)  = f(y  + n),  we  have  fn(y)  € F 
and  hence 

gjf (y  + n),y]  < gx[f (y),y] . 5.3 

Integrating  equation  5.3  with  respect  to  dm(y)  and  dm(n),  we  get 

P < P = A f gx  [f  (y),y]dm(y)  5.4 

for  all  f(0  C F. 

We  shall  assume  in  what  follows  that  the  gunner  can  reduce  noise, 
or  will  reduce  noise,  or  should  not  be  assumed  not  to  be  able  to  reduce 
the  noise  to  essentially  zero.  In  any  event,  we  can  put  an  upper  bound 
on  his  performance  by  assuming  that  he  is  in  a noiseless  situation. 

We  have  now  completely  optimized  the  problem  from  the  point  of  view 
of  the  gunner.  The  optimal  performance  from  the  gunner's  point  of  view 
is  achieved  as 

P = /g[f  Cy).y]dm(y)  5.5 


where 


g[f(y),y]  >g[f(y)»y] 


5,6 


where 


p = the  best  the  gunner  can  do 

gC*,*)  = the  conditional  probability  density  of  the  future 
position  given  the  actual  target  history.  Note 
the  change  in  notation. 

A 

f(’)  = the  optimal  fire  control  algorithm 
y = the  actual  target  history 
m(y)  = the  measure  over  the  ensemble  Y of  histories  y 

6.  OPTIMAL  TARGET  MANEUVERS 

A reasonable  course  of  action  for  the  pilot  is  to  assume  that  the 
gunner  is  smart  and  that  the  thing  for  the  pilot  (actually  for  the 
ensemble  of  pilots)  to  do  is  to  pick  g C • » • 3 and  m(y)  to  minimize  p in 
equations  5.5  and  5.6. 

Because  the  gunner  can  refire,  the  pilot  cannot  perform  only  one 
maneuver  but  must  carry  out  a series  of  maneuvers.  For  example,  if  the 
target  had  an  acceleration  limit  of  a and  initiated  an  arbitrary  but 
fixed  acceleration  after  tQ,  the  target  would  be  uniforr ly  distributed 
in  a circle  of  radius  ^jaT2  at  time  t0  + T.  A little  thought  shows  that 
this  would  be  optimal  for  a shot  fired  at  tQ,  But  if  the  gunner  were 
then  to  fire  a round  at  time  t0  + e,  he  would  be  able  to  predict  almost 
exactly  where  the  target  (now  accelerating  constantly)  would  be  at  time 
t0  + T + £. 

This  course  of  action  requires  some  sort  of  assumption  on  the 
dependence  of  g(,,0  and  m(y)  on  tQ.  Typically  one  would  assume  that 
g(<,>)  and  m(z)  were  independent  or  t0  (at  least  for  the  evasive  com- 
ponents of  the  motion).  The  impact  of  such  assumptions  on  the  character 
of  possible  maneuvers  is  discussed  in  Reference  1.  We  shall  assume  this 
stationarity  at  least  over  some  time  interval  of  the  order  of  size  of 
the  flight  time  T of  the  bullet. 


In  Appendix  B we  show  that  if 


oT2  = /{[xi-  y j (y)  ] 2 + [x2-  u2(y)]2^g(x,y)dxdm(y) 

6.1 

where 

^i (y)  = f xig (x ,y)dxdm (y) , 

6.2 

then 

* _ A 

P > 7. 

2tt  o 2 

6.3 

^S.  S.  Wolff  and  P.  R.  Schlegel,  "Fire  Control  Predictors  for 
Target  Maneuvers,"  Ballistic  Research  Laboratories  Report  No. 
April  1976.  (AD  #B011320L) 
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.... 


Further,  since  the  second  moment  is  minimized  around  the  mean,  we  have 
for  any  functions  Vj(«)  and  v2(«) 

ar2  < /{[xi-vl(y)]2  + [x2  - v2(y)]2  | g(x,y)dxdm(y)  6.4 


In  particular,  let  and  v2(>)  be  the  optimal  (r.m.s.)  linear  estima- 

tors of  the  mean  as  given  by  the  Wiener  theory.  If  a2  is  the  variance 
from  the  Wiener  theory,  then  by  equation  6.4 


and  hence 


< a2 
— rw 


A 

P > 


6.5 


6.6 


This  is  a very  useful  relationship,  a 2 may  well  not  be  easy  to 

come  by.  However,  cr2w  is  obtainable  from  just  the  correlation  function 

<KO  without  any  further  knowledge  of  the  process  (see  References  2 and  3). 

Even  though  the  optimal  linear  estimator  may  be  a poor  estimator  of  the 

mean  and/or  the  optimal  fire  control  algorithm,  we  can  still  formally 

calculate  o2  and  still  use  it  in  equation  6.6. 
rw 

If,  however,  the  pilots  assume  a Gaussian  ensemble  of  maneuvers, 
the  Wiener  estimator  becomes  the  optimal  fire  control  algorithm,  and  the 
distribution  of  prediction  errors  is 


_ [(Ax,)2*  (Ax,)2] 


no 


rw 


rw 


6.7 


where  AXj  and  Ax2  are  measured  from  the  predicted  mean.  Since  the  pre- 
dicted mean  is  the  aim  point,  and  with  our  definition  of  vulnerable  area, 
we  have  for  this  case 


P = 


A 

OTT  2 


rw 


6.8 


So  if  p is  the  best  the  pilot  could  do  in  minimizing  p,  we  have 


"Norbert  Wiener,  Extrapolation , Interpolation,  and  Smoothing  of  Stationary 
Time  Series,  The  Technology  Press  of  M.I.T.  and  John  Wiley  6 Sons,  Inc., 
New  York. 

3H.  W.  Bode  and  C.  E.  Shannon,  "A  Simplified  Derivation  of  Linear  Least 
Square  Smoothing  and  Prediction  Theory,"  Proceedings  of  the  Institute 
of  Radio  Engineers,  Vol.  38,  pp.  417-425,  April  1950. 


1 


1 . a i 

— 7~  > P > 7 

on  2 on. 


6.9 


rw 


rw 


The  left  side  obtains  since  the  pilot  can  do  that  well  by  using  a Gaussian 
process.  How  close  he  can  get  to  the  right  side  of  6.9  is  an  open  question, 
however  (see  Reference  1). 

Equation  6,9  is  especially  useful  since  we  can  incorporate  con- 
straints. For  example,  in  Reference  4 (and  outlined  in  Appendix  C)  the 
following  problem  is  solved: 

Maximize  o 2 subject  to 
rw 


E{[xl  CM)  (t)]2  ♦ [x2°°(t)]2}  = C,^2. 


The  solution  is 


rwM 


- w 


6.10 


6.11 


where  k^  is  the  largest  eigenvalue  of  a differential  equation. 

Now  this  is  as  large  as  the  pilot  can  make  o2,,  subject  to  an  r.m.s. 
constraint  on  his  Mth  derivative.  Thus 


rw 


t rtr^  „ ~ ^ ~ 2va2  „ 

rwM  rwM 


6.12 


for  all  processes  for  which  the  Mth  derivative  has  the  r.m.s,  value  of  C^., 

So  "for  any  type  of  maneuver  process  we  can  bound  within  a factor 
of  two  the  best  that  a pilot  can  do  against  the  best  a gunner  can  do  where 
the  vehicle  has  an  r.m.s.  constraint  on  its  Mth  derivative."  How  close 
the  pilot  can  come  to  achieving  the  right-hand  side  of  equation  6,12  is 
still  a matter  for  conjecture. 

7.  CONCLUSIONS 

In  addition  to  the  fact  that  equation  6.12  allows  bounds  to  be  put 
on  p,  we  can  draw  some  other  conclusions: 

• If  the  pilot  is  forced  to  use  maneuver  tactics  which  result  in 
peaks  in  g(x,z),  the  fire  control  designer  is  well  advised  to  create 
optimal  algorithms  and  design  an  accurate  and  quiet  fire  control  system 


Harry  L.  Reed,  Jr.,  "A  Simple  Model  of  Gunfire  Predictors,"  Ballistic 
Research  Laboratories  Report  No.  1775,  April  1975.  (AD  #A009507) 
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to  take  advantage  thereof.  The  size  of  the  benefits  which  might  accrue 
in  this  instance  has  not  been  discussed  in  this  paper  (except  that 
« A 

p > j ),  but  some  idea  of  possible  advantages  is  given  in  Reference  5. 

2iror 

• If  the  pilot  has  rather  complete  freedom,  he  can  achieve  the  per- 
formance specified  in  equation  6.12.  He  does  this  with  either  an  ensemble 
of  Gaussian  maneuvers  or  with  an  ensemble  with  even  a "flatter"  (if  a 
heuristic  concept  can  be  exercised)  conditional  probability  density. 

• Assuming  an  optimal  fire  controller,  there  is  a hit  probability 
associated  with  a Gaussian  target  tactic;  the  pilot  can  do  no  better 
than  reduce  this  probability  by  half. 

• A word  of  caution  is  appropriate  here.  If  the  target  is  a high- 
performance  aircraft  and  if  the  pilot  is  indeed  able  to  conduct  optimal 
maneuvers,  the  choice  of  an  optimal  fire  control  algorithm  may  be  aca- 
demic. That  is,  the  optimal  algorithm  may  be  only  the  least  poor  of  a 
set  of  poor  choices. 

• It  is  important  to  know  more  about  the  statistics  cf  targets: 

••  Will  the  maneuvers  be  mild  and  uniformly  distributed? 

••  Will  there  be  non-Gaussian  eccentricities  that  offer 
significant  gains  to  an  optimal  algorithm? 

••  Will  the  targets  maneuver  so  severely  and  unpredictably 
that  even  optimal  algorithms  would  be  inadequate  at 
certain  ranges  of  interest? 


5Harry  L.  Reed,  Jr.,  "Limitations  of  the  R.M.S.  Criterion  for  Fire 
Control,"  Ballistic  Research  Laboratories  Report  No.  1805,  July  1975. 
(AD  #A014986) 
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APPENDIX  A 


Herein  we  shall  consider  a simple  bimodal  model  that  shows  the  prob- 
lems one  can  get  into  by  using  the  mean  as  a predictor.  The  model  seems 
to  be  almost  too  simple  to  be  real,  but  there  are  at  least  two  real  ana- 
logues: a target  that  can  maneuver  hard  left  or  hard  right  immediately 

after  a bullet  is  launched,  and  two  targets  which  are  causing  a radar 
tracking  system  to  "glint”  between  them. 

Anyway,  in  equation  2.1  let 

g(x)  = j [fi(x  - 1)  + 6(x  + 1)] 

v (a  - x)  = 1 |a  - x | < e 

= 0 otherwise 

b(0  = «(•) 

We  note  that  this  g is  independent  of  z,  so  we  can  take  f(«)  to  be 
independent  of  z.  Integrating  and  changing  some  dummy  variables,  we 
have 


f - 1 ♦ c 

-f  ♦ 1+  e 

j (c)dc  ♦ J 
f - 1 - e J 

' f + 1 - e 

j(c)dc 


and  if  0 < e < 1,  the  regions  of  integration  do  not  overlap,  so 


hfj 


(c)dc  = . 


So  for  any  prediction  algorithm  and  any  dither,  p is  bounded  by  y . 

An  optimum  algorithm  would  be  one  for  which  p equals  y.  A 
class  of  optimum  algorithms  is: 

Let  f = 1 with  probability  a 

= -1  with  probability  1 - a. 

Then  we  would  have 


p ■ 


j (c)dc  ♦ 


r2  + e 

^ / j 

J 2-e 


j(c)dc 


♦ (1  - «) 


1 

2 


e 

j (c)dc 
e 


We  can  take 

j(0  = 6(0 

and  have 

P ■ 7 « ♦ T 0 ' °>  ' 7 • 


Note  that  this  can  be  a completely  deterministic  algorithm  by  choosing 
o * 0 or  o = 1.  Clearly  as  e goes  down  in  size  j (c)  must  concentrate 
its  probability  mass  at  0,  2,  or  -2  if  we  are  to  achieve  p = 1/2. 

If,  on  the  other  hand,  we  were  to  use  the  mean  as  the  predictor, 
we  would  set  f ■ 0 and  have 


P 


1 

2 


/ 


-1  ♦ e 

j (c)dc  + 

-1  - e 


1 

2 


j (c)dc. 


Now  we  could  use  delta  functions  at  1 and  -1  to  define  j(0»  but  this 
would  just  be  another  way  of  getting  an  optimum  prediction  algorithm. 
More  usually,  j(0  would  be  characterized  by  something  that  at  least 
would  be  unimodal  and  symmetric,  in  which  case  we  would  have 


Now 


and 


so 


1 
/. 


j (c)dc  < j(l  - e)2e 


e 

1 - e 

0 

1 -e 


j (c)dc  < j 


J(c)dc  > (l-e)j(l-e) 


2c  e 

P - 2(1  - t ) " 1 - e 
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Msaamm 


I 

Therefore  p is  small  when  e is  small  rathei  than  i or  any  value  of  c 
as  is  the  case  for  the  optimum  predictor. 

Nevertheless,  dither  helped  the  mean  estimator  since  we  would  have 

P ■ 0 

for  f 0 and  j ( • ) = 6 (.  • ) . 


v r 

prscsdiic  PiaiA'<UMc^wr  filmed 


APPENDIX  B 


In  this  section  we  shall  derive  equation  6.3.  For  convenience  we 
shall  just  consider  the  two-dimensional  probability  distribution  function 
g(x)  where  we  have  suppressed  the  dependence  on  y.  We  shall  also  con- 
sider for  convenience  that  g(x)  has  zero  mean. 


or2  = /(Xl2  ♦ x22)g(x)dx1dx2. 


We  can  change  to  polar  coordinates  and  have 


} = J r2g(r,6)rdrde 

1 * J g(r,0)rdrd6. 


and  we  can  define 


g\r) 


* 2 IT 

= I g(r,e)de. 

J 


With  this  definition  we  have 


Now  let 


and  let 


J g*(r)rdr  = 1 

>r2  - /’rV 


(r)dr 


6 = sup  g(x). 


2 1 

Y = 


0 < g*(r)  < 2ttB, 


We  can  write 


(r)dr 


5r2  ■ J r3g*(r)dr  ♦ J'  r3g*(r)dr 

^ f r3g*(r)dr  ♦ y2  f rg*(r)dr 
Jo  ' Y 

* J r3g*(r)dr  ♦ Y2  | 1 - rg*(r)drj 

■ f r3g*(r)dr  ♦ y2  j J -~z  rdr  - J'  rg* 

■ f r3g*(r)dr  * y2  f [2*6  - g*(r)]rdr 

^ f r3g*(r)dr  ♦ f [2wB  - g*(r)]r3dr 
^ 0 J O 


(r)dr 


This  last  step  uses  the  fact  that 


[2tr6  - g*(r)]  ^ 0, 


Finally  we  have 


•S  r\ 


2ir6r3dr 


. ier_  _ i 


Now  we  turn  our  attention  to 


where 


P * A J g(y)dm(y) 
aT2  = J ar2(y)dm(y) 
or2(yi  = J I x - u(y)  |2g(x,y)dx 


g(y)  * sup  g (x , y) 
x 
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From  what  we  just  proved  we  can  write 


2>  (— 
J 2irg  I 


dm(y). 


2*g  (y) 


Now  we  have 


a P Sr 
r v 


2irg  Cy) 


dm(y)  / VA8(y)  | *“(y) 


which  gives 


- f l/  ■ 1 — VAg  (y)  dm(y) 
./  r 21rg(y) 


by  a version  of  Schwartz's  lemma  and  the  fact  that  m(y)  is  a probability 
measure.  This,  of  course,  gives  the  final  result 


P * 


that  appears  in  equation  6,3, 


If  we  had  been  considering  a one-dimensional  situation,  we  would 
have  had 

sup  g(x)  >■  


2v/3*  a 

x 

where  L is  a "vulnerable  length."  We  would  also  have  as  the  analogues 
of  equations  6.7  and  6.8 


(Ax)2 


. 2 a2 

i e xw 


.iM 
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PRSCEDItG  PA3 ii&ANK-NOT  FILMED 


The  argument  in  Reference  4 that  is  referred  to  in  Section  6 of 
this  paper  goes  as  follows. 

Consider  the  one -dimensional  problem.  The  optimal  (r.m.s.)  linear 
predictor  is  given  by 


x(t^T)»^]  x^(t)  “T  + /*  (t  - t)f  (£)d£ 

m*0  0 

where 

E[xW(tj)xW(t2)]  =*M(ti-t2) 

C-l 

C-2 

and 

V°>  - Si2  < - • 

C-3 

Now  we  can  write 

VT)  = f 

Jo 

C-4 

ik|(t)  = 0 for  t < 0 

C-5 

S <2-f  V‘i!dt' 

Jo 

C-6 

If  we  minimize  the  expected  value  of  the  square  error  of  equation  C-l 
by  picking  f ( • ) , we  have  for  the  residual  error 

r T r r S M_j  "J  2 

■ / ds  / otHIt  Vr>dr  c-7 

*o  L*/0 

The  pilot  of  the  target  can  maximize  e2  by  choosing  optimal  statis- 
tics represented  by  subject  to  C-6.  This  maximization  gives 

£2  * w 

C-8 

where  k.  is  the  largest  eigenvalue  of  the  differential  equation 
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*»>> 

= 0 

L = 0 to  M - 1 

C-9 

Cm 

= 0 

L = 0 to  M - 1 

C-10 

vr* 

= (- 

C-ll 

Likewise  in  two  dimensions,  if 

■ eJ^’w  2 * *2(M)(t)]2j 

then 

°rwM  = Si  Si  - 
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